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Abstract 



We present two examples of how single-molecule experimental techniques applied to biological systems can give insight 
into problems within the scope of equilibrium and nonequilibrium mesoscopic thermodynamics. The first example is the 
mapping of the free energy landscape of a macromolecule, the second the experimental verification of Crooks' fluctuation 
theorem. In both cases the experimental setup comprises optical tweezers and DNA molecules. 
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1. Introduction 

Mesoscopic thermodynamics (also known as thermody- 
namics of small systems fV) deals with the intermediate 
scale between the microscopic and macroscopic level. At 
such scale the typical number of microscopic components 
is much larger than 1 but much smaller than Avogadro's 
number, the energy exchanged between the system and its 
environment is of the order of few k^T, and fluctuations 
of thermodynamic quantities about their expected values 
are both large and meaningful, that is, they are easily de- 
tected and carry information about the structure of the 
system itself. 

Many of the systems investigated in biophysics are small 
I in the sense defined above. This is why the recent advances 
in single-molecule experimental techniques that allow 
the manipulation of biological nanostructures, such as pro- 
teins or nucleic acids, provide with a valuable investigation 
tool not only the biologist, but also the physicist interested 
into the mesoscale world. 

The goal of this paper is to offer a taste of the po- 
tentiality of these techniques; it is addressed primarily 
to physicists not already familiar with the subject. We 
have selected, among many active research fields, two top- 
ics we have been working on using optical tweezers in our 
laboratory in Barcelona. The former belongs to equilib- 
rium thermodynamics: in Sec. [3] we discuss how to map 
the free energy landscape of a biological macromolecule, 
specifically a long DNA hairpin. The latter is a problem 
of nonequilibrium thermodynamics: the experimental ver- 
ification of Crooks' fluctuation theorem, to be treated in 
Sec. [4] Before delving into the applications, however, we 
sketch in the next section a brief outline of the experimen- 
tal setting. 
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2. Optical tweezers 

Optical tweezers [3] use the pressure of light radiation to 
measure and apply force to small transparent objects, such 
as microscopic spherical particles (usually called beads) 
made of polystyrene, latex or silica. A bead illuminated 
by a laser beam undergoes two forces: one is proportional 
to the gradient of the intensity of light, the other is the 
scattering force due to the light reflected on the bead sur- 
face. If the gradient and scattering force are equilibrated, 
an optical trap is formed. This can be achieved by using a 
well focused Gaussian beam with high numerical aperture. 
The measurement of the deflection of the laser beam due 
to its interaction with the bead makes possible the evalua- 
tion of the applied force, once the principle of conservation 
of light momentum is taken into account (see Fig.[T|L). 

Our apparatus (see Fig. [T|d) is a dual counterpropa- 
gating beam miniaturized optical tweezers |4], with fiber- 
coupled diode-lasers that produce a piezo controlled mov- 
able optical trap. It is possible to apply and measure forces 
in the range 1 100 pN. The instrument also measures 
variations of the position of the optical trap's center. The 
resolution is 0.1 pN in force and 0.5 nm in position, and the 
sampling rate is 1 kHz. The experiments are performed in 
a fluidic chamber so that the molecule of interest is sur- 
rounded by water buffer. A thoroughly detailed descrip- 
tion of our experimental apparatus is due to appear in a 
forthcoming paper [5 , while an effective theoretical model 
for the characterization of the bead in the optical trap and 
the handles that hold the molecule can be found in Ref . [B] . 

The DNA molecule is labeled with biotin and digox- 
igenin groups in both ends. These labels stick to anti- 
gen coated polystyrene beads (sized few micrometers). A 
pulling experiment is set when one bead is in the trap 
and the other is held in a fixed position by a micropipette 
(see Fig. A canonical experiment consists of repeated 
pulling and relaxing cycles during which the two beads are 
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Figure 1: Optical trapping, a In the upper picture, a transparent 
bead is trapped by a focused laser beam. When the bead is cen- 
tered, the incoming laser beam has the same linear momentum as 
the outcoming one. The bead feels no force. In the lower picture, an 
off-center bead undergoes a restoring force. The laser beam is de- 
flected and the bead feels a force so that the total linear momentum 
is conserved, b Sketch of the instrument. The instrument consists on 
two symmetrical optical paths that produce one single optical trap. 
Laser beams are produced in laser diode moduli (a). The light is 
conducted to the so-called wigglers (c) through an optical fiber (b). 
The wigglers redirect the light using piezo actuators that change the 
position of the optical trap. A small amount of light is split from 
the main beam and collected by a Position Sensitive Detector (PSD) 
that records the position of the optical trap (d). The remaining beam 
is used to form the optical trap in the fluidic chamber (f) by means 
of a microscope objective (e). The outcoming beam is collected by 
the other objective and evaluated by another PSD that measures the 
force (g). The optical trap and the whole experiments are observed 
using a CCD camera (h). c Experimental setup. Each strand of 
the DNA molecule to be unzipped is bonded to equal length dsDNA 
handles and these are tethered to the bead. The upper bead follows 
the optical trap when it is moved upwards and the force transmitted 
along the handles produces the unzipping of DNA [T] |8] |9] I10| 
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Figure 2: Nearest-Neighbor model, a The formation energy of the 
molecule depends on the formation energy of each single base-pair. 
Since there are 4 types of bases (Adenine, Guanine, Cytosine and 
Thymine), 16 different NN base-pair interactions are possible. How- 
ever, due to symmetry reasons, they get reduced down to the 10 
independent parameters highlighted in the table, b The description 
of the experiment is broken down into different elements. Each ele- 
ment contributes to the total distance of the system according to its 
extension at a given force. The total energy of the system is the sum 
of the contributions of the separate elements. 



periodically separated and brought back together. The ra- 
pidity of the process is characterized either as puUing speed 
(expressed in nm/s) or as loading rate (measured in pN/s). 
Data are collected in the form of force-distance curves (see, 
for example, Fig. [6| , from which a wealth of information 
about the molecule can be mined. 



3. Mapping the free energy landscape of DNA 
molecules 

Using the apparatus described in the previous section, 
we have performed mechanical unzipping experiments on 
long (thousands of base-pairs) DNA molecules [5]. These 
experiments are well described by the Nearest-Neighbor 
(NN) model for nucleic acids [11 , which predicts the free 
energy of formation of a dsDNA (double-stranded DNA) 
from two complementary sequences of ssDNA (single- 
stranded DNA). According to the NN model, the free en- 
ergy of formation of the j-th base-pair ej depends not only 
on the bases that appear in position j, but also on the 
nearest-neighbor base-pair in position j + 1 (see Fig. [2^). 
The total free energy of formation of a molecule is simply 
given by the sum of all individual contributionsj^ 
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(1) 



where the discrete variables aj take values in the set of 4 
canonical couplings between the bases: 

CTj e {AT, TA, CG, GC} . 

Such free energy is explicitly sequence-dependent. 

The additive nature of the free energy (jlj allows us to 
calculate the free energy G^sW) of intermediate states, i.e., 
when the first n base-pair bonds are broken. An interesting 
application of this model is that a DNA sequence can be 
specifically designed in order to obtain the desired free 
energy profile. 

Provided that the process is in quasi-equilibrium and the 
force / along the system is homogenous at any moment, 
the total free energy is simply the sum of the contribu- 
tions from all the components: the trap, the handles, the 
double- and single-stranded DNA. The energy i?b of the 
bead in the trap is approximately quadratic in the distance 
a;b of the bead's center from the trap's focus: 



-ky,X 



2fcb 



(2) 



where k\, is the stiffness of the trap, function of the power 
of the laser and the properties of the bead (shape, size, and 
material), and the latter equality follows from / — khX^. 



^Actually, one needs to distinguish a loop contribution at the end 
of the chain. 
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Open basepairs 

Figure 3: Free energy landscape of a 2252 base-pair (bp) molecule at 
fixed xtot- a Several free energy vs. number of open base-pairs curves 
are plotted for different total distances. At a fixed total distance, 
the energy landscape shows a minimum of energy which is assumed 
to be the most probable state. The inset shows a detailed view 
where the roughness of the free energy landscape can be appreciated. 
The global minimum is surrounded by other local minima that may 
coexist with it. 

The elastic free energies of the polymers are calculated by 
integration of the appropriate force vs. extension curves: 

Gi(xi)= / ' Fi{x)dx (3) 

JQ 

where i is either s for the ssDNA (modeled by a Stretched 
Freely- Jointed-Chain [12]) or h for the handles (modeled 
by a Worm-Like-Chain [TB]V 

The explicit calculation depends on the statistical en- 
semble, which is determined by the control parameter of 
the unzipping experiments. For optical tweezers, the con- 
trol parameter can be either the total extension xtoti equal 
to Xb + 2a;h + 2xs, see Fig. ^p) or the applied forccjj/. In 
the former case the free energy is: 

G{xtot,n) = E],{xb) + 2Gh(xh) + 2Gs(a;„ n) + Gds(?i) • (4) 

As for the constant force ensemble, the calculation of 
the total free energy H(f, n) comprises the same elements 
that appear in Eq. ([I]): 

H{f, n) = Sb(/) + 2i7h(/) + 2F,(/, n) + GM , (5) 



^Magnetic tweezers are ideal setups to control force 1101 . Also 
specifically designed tweezers setups with zero stiffness regions |14| 
can operate at controlled force. Force-feedback systems are not ideal 
constant force systems as they introduce other sort of noise effects 
due to the limited feedback frequency (typically around 1 kHz). 
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Figure 4: Experimental data of hopping transitions between 3 basins 
that coexist. Although the effects of some intermediate states are 
detected, they can not be identified because their transitions are 
masked by thermal fiuctuations. The theoretical calculation of the 
free energy landscape allows us to recognize 3 minima or basins that 
are separated by barriers that are in the range 5-^10 k^T, as shown 
in the inset. The free energy difference among them is small enough 
to experimentally observe the coexistence. 

where is written in terms of the applied force and the 
free energy of the polymers is given by the Legendre trans- 
form of Eq. (|3|: 

mf) = - X.indf , (6) 
Jo 

where i stands for either s or h and Xi(f) is the inverse of 

Equations (|4| and (|5| are the basis for the study of the 
equilibrium thermodynamics of DNA unzipping. Several 
important properties can be investigated; here we focus on 
the reconstruction of free energy landscapes. 

Figure [3] shows a plot of the numerically computed [15j 
free energy landscape at fixed total distance. The tun- 
ing of the control parameter xtot changes the shape of the 
landscape, and the stability of the states characterized by 
n is being modified. In a typical unzipping experiment, 
where xtot is slowly increased at a fixed rate, the DNA 
molecule passes through a succession of states character- 
ized by the number of broken bonds n. The unzipping is 
not a correlative process in which each base-pair bond is 
broken before the next one. Instead, it is a cooperative 
phenomenon in which the base-pair bonds are broken in 
groups. This agrees with a rough landscape where min- 
ima in the free energy abruptly change from state rii to n2 
(where n2 — rii is usually larger than 1) as Xtot increases. 

The experiments also show coexistence of states in which 
hopping between them can be observed. The residence 
time of the states is about 1 s, so the time resolution of 
the instrument (1 ms) allows us to identify this process. 



3 




Figure 5: Free energy landscape of a 2252 base-pair (bp) molecule at 
fixed /. Three free energy vs. number of open base-pairs curves are 
plotted for different forces which induce different tilts of the land- 
scape. There is a coexistence force (/ = 16.5 pN) at which the mean 
tilt is zero. In this situation, many states with different number of 
open base pairs coexist. 



This phenomenon is well described by the roughness of the 
free energy landscape, too. Figure |4] shows a magnified 
region of the free energy landscape when the coexistence 
of 3 states is observed for a finely tuned value of a;tot- 

As can be seen in Fig. |5] similar but slightly different 
aspects and details are emphasized by the free energy land- 
scape of the system at fixed force. In this case, the force 
induces a tilt in the free energy landscape. From this de- 
scription, we can deduce that there is a coexistence force 
that induces no tilt in the force landscape. Under this sit- 
uation, there are a lot of minima along the molecule that 
coexist and may be observed. If the force is increased, the 
tilt is induced, the states with higher number of broken 
bonds become more stable and the molecule begins unzip- 
ping. On the contrary, if the force is reduced the molecule 
is re-zipped. Interestingly, the tilt can be adjusted by tun- 
ing the constant force applied on the molecule. 

The features of the system we are exploring allow us to 
design experiments to study the diffusion in an unidimen- 
sional rough landscape. A particular sequence of DNA 
can also be synthesized in order to produce desired peaks 
and valleys. The combination of the experimental mea- 
surements and the whole theoretical description has been 
fruitful and allows us to explore the equilibrium and close- 
to-equilibrium properties of DNA base-pair interactions. 

4. Experimental test of Crooks' fluctuation theo- 



The understanding of nonequilibrium thermodynamics 
of small systems [ W\ experienced in the last ten years a 
season of rapid progress, whose milestones are the theo- 
retical results collectively known as fluctuation theorems 
|17j and their experimental verifications. Here we are con- 
cerned with one such theorem, due to Crooks 



Figure 6: Some representative force-distance curves. Here the "fast" 
cycle of unfolding/refolding is carried on at 400 nm/s (corresponding 
to a loading rate of 18.4 pN/s), and the "slow" one at 40 nm/s (or 
1.84 pN/s). 



Let us consider a small system immersed in a thermic 
bath at temperature T (for instance, our setup of optical 
trap -I- beads -I- handles + DNA molecule). Let us say 
we can manipulate the system by varying the value of a 
control parameter A (for us, the total distance Xtot between 
the trap and the pipette). Other parameters that identify 
the state of the system and cannot be affected directly, we 
group them together in a variable s(A). In our example, 
s may be a discrete variable assuming value if the DNA 
hairpin is closed and 1 if it is operj^ If the system is in 
thermodynamic equilibrium, then the probability density 
of s follows the Maxwell-Boltzmann distribution. Now we 
choose an experimental protocol Xpit), which we label as 
'forward', that starts in equilibrium at time t ~ with 
A = Aq and ends out of equilibrium at time t = t with 
A = At-. Another protocol, labeled 'reverse', also starts 
from equilibrium and is such that Ar(<) = Af(t — t). 

Let the work W be defined as the energy that we feed 
into the system (the whole experimental system, trap in- 
cluded) throughout the nonequilibrium process of varying 
A. We indicate with VpiW) and 7^R,(Vt^) the work proba- 
bility densities along the forward and reverse process, re- 
spectively. Crook's theorem states that, provided that the 
microscopic dynamics always obeys the detailed balance 
condition iTOl. 



Vf{W) 
Vn{-W) 



'W-AG\ 



(7) 



where fee is the Boltzmann constant, and AG is the equi- 
librium free energy difference 



AG = G(A„s(Ar))-G(Ao,s(Ao)). 



(8) 



■^Here we are considering a DNA molecule short enough that it 
opens all at once, without the intermediate states considered in sec- 
tion |3] 




Figure 7: Work distributions 7'f(W) (unfolding) and 'Pr(— W) (re- 
folding). The vertical lines highlight the crossing points. The dis- 
tance between these two estimates of the free energy is about 1.1 
k^T. The distributions are evaluated from 561 fast and 228 slow 
unfolding/refolding cycles. 



From Eq. with elementary manipulations, we get 

exp(-/3AG)7'R(-T4^) = Vf{W) exp(-/3iy) , (9) 

where p is as usual the inverse of k^T. Now we can mul- 
tiply both sides by a generic function (j){W) and integrate 
over W. The result can be arranged in the form 



exp(-/?AG) = 



(10) 



where the angular brackets (• • • )f(r) stand for an average 
over all possible realizations of the forward (reverse) pro- 
tocol. The simplest possible choice, 4>{W) = 1, yields the 
well-known Jarzynski identity |20| . To the purpose of mea- 
suring free energies, however, there is a more convenient 
option. 

In practice, one estimates the averages that appear in 



Eq. (10 1 using a finite number of experimental events. 



Equation ( 10 1 can therefore be interpreted as the defini- 



tion of an estimator of AG. It turns out that the statistical 
variance of such estimator is minimized by Bennett's func- 
tion ED 



1 



exp[P{W - AG)] 



(11) 



where n-p, tir is the number of forward or reverse events, 
respectively. 

Following the steps of the first experimental test of 
Crooks' theorem [22j, we have performed pulling exper- 
iments with optical tweezers on a short (20 bp) DNA 
hairpin, characterized by a two-state behaviour [23 . In 
general, when one is interested in studying some physical 
property of a particular molecule, it is advisable to collect 
data from as many specimens as possible. Here, however, 
our goal is to illustrate the validity of Eq. and the 



Figure 8: Bennett's so-called "acceptance ratio" method. The nearly 
horizontal lines are the functions z^{u) — zp{u) for fast and slow 
trajectories. They intersect the line y = u at 413.53(20) k^T and 
412.91(14) fceT, respectively. 



sake of clarity would be poorly served by the variability 
brought about by the inevitable individual differences in 
a large sample. We use therefore data taken at two dif- 
ferent pulling speed (fast, 400 nm/s, and slow, 40 nm/s), 
but from the same specimen. In this way we are assured 
that we are always working with the same object, and can 
concentrate on the issue of how precisely the Crooks fluc- 
tuation relation can be checked in actual experiments. 

In Fig. [6] we show few typical force-distance curves /(A). 
The 1 pN jump in the force is the signal of the unfolding or 
refolding transition. It may be observed how the fast cycle 
exhibit larger hysteresis, while in the slow process it is not 
infrequent to have multiple transitions between the closed 
and the open states. As the pulling speed is reduced, the 
so-called transient violations of the second law [24; . that 
is, cycles in which the energy expense is negative, also 
become more probable. 

From the force-distance curves we can compute the me- 
chanical work performed on the system 



W : 



/(A)dA. 



(12) 



Figure [7] shows the estimated work distributions thus ob- 
tained. The work probability density functions have been 
evaluated according to a histogram-free method recently 
proposed in Ref. one nice advantage over tradi- 

tional, histogram-based methods is that the theoretically 
expected asymmetry of the work distribution [26 is evi- 
dent, at least in the "fast" case. 

Equation ([7| implies that the graphs of Vf{W) and 
VrI—W) cross each other for W = AG: the work dis- 
tributions themselves are dependent on the pulling speed, 
but the crossing point is not. In our example, the dis- 
tance between the crossing point of the fast distributions 
is about 1.1 k^T ~ 0.6 kcal/mol larger than the crossing 
point of the slow ones. However, we should observe that 
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research. In this paper we have illustrated how they 
can be equally useful for the physicist interested in the 
mesoscale phenomenology. Easily synthesized biological 
macromolecules are an ideal laboratory to verify our com- 
prehension of equilibrium and nonequilibrium thermody- 
namics of small systems. 
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Figure 9: Plot of \n{Vv (W)/7'r(- VF)) as a function of W. The least 
squares linear fit slope is 0.96(7) for the "fast" data and 1.0(1) for 
the "slow" ones. 



the errors on the "slow" distributions are large, besides, 
this is hardly the best way of estimating AG. As explained 
earlier, the most convenient way to extract a free energy 
estimate from Eq. ([t]) is to use Bennett's choice (111 into 
Eq. ([T0|. 



In practice, we define the functions 



zf(u) ~ In 
zr(u) = In 



exp(-/3W^) 
1 + exp[/3(M^ - u)\ I p 
1 

1 + exp[-/3(M^ + u)\ 



(13a) 
(13b) 



where we are taking advantage of the fact that for us np = 
nR, and plot the function, that should be approximately 
constant, zr(u) — zf(u) (seeFig.|8|. The intersection with 
the line y — u \s the optimal estimate of AG. To this 
number we can attach a standard deviation, as explained 
in Ref. [SZ]: in our example, we find AG = 413.53(20) 
k^T from the "fast" data and AG = 412.91(14) /cbT from 
the "slow" ones. The discrepancy between the two values 
has been reduced to about 0.6 k^T ~ 0.4 kcal/mol, and 
is consistent with the hypothesis that they are in fact two 
estimates of the same quantity, as required by Crooks' 
theorem. 

Another way of testing Eq. ([t]) is to directly plot the 
quantity ln(pF(l^)/^R(-T^)), as in Fig.|9l If the Crooks 
fluctuation relation is satisfied, then the curves must be 
straight lines with slope 1 which cross the VF-axis at /3AG. 
The slopes of the linear fit that we find in this way, 0.96(7) 
and 1.0(1) for the fast and slow case, respectively, are fully 
consistent with the theoretical expectations. 
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5. Conclusion 

Single-molecule experimental techniques such as the 
atomic force microscope and optical or magnetic tweez- 
ers are now established as a routine tool in biophysical 
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